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Abstract 

Multipole electron modes beyond the Mie plasmon in atomic clusters are investigated within 
the time-dependent local density approximation theory (TD-LDA). We consider the origin of the 
modes, their connection with basic cluster properties and possible routes of experimental obser- 
vation. Particular attention is paid to infrared magnetic orbital modes, scissors and twist, and 
electric quadrupole mode. The scissors and twist modes determine orbital magnetism of clusters 
while the electric quadrupole mode provides direct access to the single electron spectra of the 
cluster. We examine two-photon processes (Raman scattering, stimulated emission pumping and 
stimulated adiabatic Raman passage) as the most promising tools for experimental investigation 
of the modes. 
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I. INTRODUCTION 



Besides the dominant dipole Mie plasmon, many other electron modes, both electric and 
magnetic, can exist in atomic clusters These modes are known in diverse many-body 
systems (atomic nuclei, quantum dots, dilute gas of trapped fermionic and bosonic atoms, 
etc) . It is a demanding problem to observe them in atomic clusters. The multipole electron 
modes represent an essential part of electron dynamics in both collective and electron-hole 
domains. They are connected with basic cluster properties and thus can serve as an effective 
tool for their investigation. 

We will pay the main attention to three infrared modes: scissors magnetic dipole (Ml), 
twist magnetic quadrupole (M2) and electric quadrupole (E2). These modes are connected 
basic properties of atomic clusters. Besides, they are the strongest spin-saturated modes 
beyond the Mie plasmon. 

The scissors and twist modes are of orbital magnetic character 0, [J U, . In general, the 
scissors mode exists only in deformed systems. It dominates the Van Vleck paramagnetism 
3, 0| and can result in dia-para magnetic anisotropy in particular light atomic clusters 
6(. The twist mode is the strongest magnetic orbital mode in spherical clusters where the 
scissors mode vanishes 5(. It is mainly generated by transitions between electron levels 
with maximal orbital moments. Altogether, the scissors and twist modes are fundamental 
sources of orbital magnetism in spin-saturated clusters. In deformed clusters, the scissors 
Ml is coupled to electric quadrupole E2 mode. 

The infrared E2 mode is most interesting in two particular cases of deformed clusters: 
free light deformed clusters and embedded oriented rods (strongly prolate large clusters). In 
the first case, the E2 mode is reduced to a few electron-hole excitations driven by cluster 
deformation. As was recently demonstrated, these excitations allow to determine the mean 
field spectra of light clusters pj. Being sensitive to cluster structure, these spectra can deliver 
important information on diverse cluster features. And last, but not least, the infrared 
electron-hole quadrupole excitations have a good chance to be observed in particular two- 
photon processes. 

Infrared E2 modes in oriented silver rods embedded into glass matrices (IJ 0] represent 
another useful example of non-dipole electron motion in clusters. These modes are collective 
and in principle can be studied in Raman scattering (experiments of this kind are in progress 
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[icl]). Because of the extreme axis ratio, rods represent a unique sample for investigation of 
both E2 and Ml (scissors) collective motion. 

Being non-dipole, the infrared multipole states cannot be populated by one-photon tran- 
sitions which excite exclusively dipole modes. Thus one has to use two-photon processes 
(TPP) where the target state is populated via an intermediate dipole state by two (absorp- 
tion and emission) dipole transitions. The dipole plasmon or isolated dipole states can serve 
as the intermediate state. Some of the processes, Raman scattering (RS), stimulated emis- 
sion pumping (SEP), and stimulated adiabatic Raman passage (STIRAP), are widely used 
in atomic and molecular spectroscopy to populate non-dipole modes but their applications 
to clusters are still very limited and, at a first glance, even questionable. The problems are 
caused by particular clusters properties (dense spectra, broad level structures, harmonicity 
of collective modes, short lifetimes, non-radiative decay channels, etc) which can hamper 
TPP. Nevertheless, we will show that TPP listed above can be applied to clusters and deliver 
valuable information on cluster properties (see preliminary discussion and first estimations 
in [3]). Specific requirements to TPP applications will be considered in detail. 

The manuscript provides a survey of explorations from a theoretical perspective. To 



that end, we employ the microscopic time-dependent 



theories in the linear regime 
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ocal-density-approximation (TDLDA) 



131 ] . The main aim of our study is to 



llj and beyond it 

outline non-dipole electron modes in clusters and encourage their investigation in TPP 
experiments. 

The paper is outlined as follows. In Sec. |Hj we give a brief overview of infrared multipole 
modes, clarify their origin and connection with basic cluster features. In Sec. IIII| the two- 
photon processes RS, SEP and STIRAP are inspected as possible routes of experimental 
study of the modes. Mainly free light deformed clusters are considered. A summary is given 
in Sec. llVl 



II. HIERARCHY OF MULTIPOLE MODES 
A. Classification scheme 

Multipole modes are oscillations of multipolarity A/i which are typical for finite many- 
body systems. It is useful to distinguish collective and non-collective multipole modes. In 
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the first case, one deals with a superposition of many elementary excitations which form 
a coherent motion while in the second case the mode represents a mixture of only a few 
elementary excitations. In mesoscopic systems with 10 — 10 4 particles, we have an interme- 
diate (and maybe the most complicated) situation when a large number of small elementary 
components contribute coherently and thus form the collective part of the mode while a 
few dominant components determine its particle-hole structure. Just such systems will be 
considered in the present study. More specifically, we will deal with metal clusters where 
valence electrons move freely in a common mean field. The single-particle levels of the mean 
field are bunched into quantum shells [3, [l^ which constitute a key point in understanding 
the nature of multipole modes and their classification. 

A useful sorting scheme for single-particle levels of valence electrons in metal clusters is 
provided by the three-dimensional harmonic oscillator The levels of a spher- 

ical cluster are sorted in perfectly degenerate bunches (major quantum shells) which are 
characterized by the principle quantum number M = 0,1,2,.... The shells are separated 
by appreciable energy gaps and every shell involves only states of the same space parity 
7r = (— 1) . This oscillator picture is well fulfilled in light clusters and provides still a good 
approximation in medium and heavy ones. In axially deformed clusters, the single-electron 
levels are characterized by Nilsson-Clemenger quantum numbers v = [\fn z A] where the 
principle shell number M = n z + 2n r + A is expressed through the numbers of nodes in 
radial (n r ) and symmetry axis (n z ) directions and projection A of the orbital moment onto 
the symmetry axis [lfil ]. 

Following this scheme, excitations of valence electrons are characterized by the AjV value, 
the difference in shell for the dominant leh jumps. Quadrupole excitations have even parity 
and are collected into the branches AJV = and 2. The excitations with larger AA/" are 
weak and can be neglected. The branch AAf = has low excitation energy. It exists only 
in deformed systems (with partly occupied valence shell) and vanishes in spherical systems 
(with fully occupied valence shell). In the next subsections, we will use this sorting scheme 
for explanation of the origin and properties of quadrupole E2, scissors Ml and twist M2 
modes. 
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FIG. 1: Photoabsorption dipole cross section in light axially deformed Na clusters. Quadrupole 
and hexadecapole deformation parameters are indicated in boxes. RPA results are given as vertical 
bars (in eV A 2 ) and as a strength function smoothed by the Lorentz weight with the averaging 
parameter 0.25 eV. Separate contributions to the strength function from the \(i =10 and 11 dipole 
branches (the latter is twice stronger) are shown by dash curves. The experimental data (triangles) 
from j^lj are given for the comparison. 

B. Calculation scheme 

The electron cloud is described by density-functional theory at the level of the local- 
density approximation (LDA) for the ground state and time-dependent LDA (TDLDA) for 
the excitations, using actually the functional of [l^. The ionic background of the cluster is 
approximated by the soft jellium model allowing for quadrupole and hexadecapole deforma- 
tions Q] . In the analysis given in this section, the infrared modes stay in the regime of 
small amplitudes. We thus employ the linearized TDLDA, often called the random-phase- 
approximation (RPA). The actual implementation for axially symmetric clusters is explained 



)prox 

u 



in |TJj . The reliability of the method has been checked in diverse studies of the Mie plasmon 
in spherical |3] and deformed jll[ I^J clusters. As an example, Fig. [T] demonstrates quite 
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good agreement of our results with the experimental data [21] for the dipole plasmon in light 
deformed clusters at room temperature. 

The proper choice of the light clusters to be studied is very important, i) The clusters 
should be small enough to possess a dilute and non-collective infrared spectrum. Only then 
the spectrum can be resolved and unambiguously related to the single-particle levels, ii) 
Since infrared modes are mainly induced by cluster deformation (see discussion in the next 
subsection), the clusters with a strong deformation (either prolate or oblate) are desirable. 
The simplest case of axial shape is most suitable for the analysis, iii) Shape isomers exhibit 



different single-electron spectra [Hi, 120, |22j , which can smear out the low-energy spectral 



lines. The heavier clusters, the more isomers 
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. So, light clusters with one dominant 



equilibrium shape are preferable. Between them, we should choose the clusters whose ground 
state and first isomers have the similar (prolate or oblate) shape. Thus we will minimize 
the spectral blurring, iv) The jellium approximation is certainly rough for description of 
details in light clusters, which establishes a lower limit for the cluster size. The results for 
the smallest samples may not reach a quantitative level, but they are still useful for the 
first consideration. As is demonstrated in Fig. ^ even in the lightest clusters the jellium 
TDLDA sufficiently well reproduces the basic characteristics of the Mie plasmon (average 
energy, principle gross- structure, magnitude of the deformation splitting of the resonance), 
see for discussion [?} 0] . Such accuracy suffices for our present survey. For this reason, as 
well as for the sake of simplicity, our study will exploit the jellium approximation. 



C. Quadrupole mode 

We start our analysis with quadrupole excitations in light deformed free sodium clusters. 
The hierarchy of E2 modes is illustrated in Fig. 121 in terms of normalized quadrupole 
photoabsorption [ill cr(E2,u;)/N c ~ Z^OI^^/JO) 2 ^^^ — uj) calculated within RPA (here 
ujj and is the energy of RPA j-state and 77 is the Lorentz weight function). Though E2 
modes are not observed in the photoabsorption, it is instructive for a first overview. Let us 
first consider the strong quadrupole resonance appearing at high frequencies in the range 
2-4 eV. In heavy clusters this resonance is associated with the quadrupole plasmon. The 
resonance is mainly formed by E2 transitions over two major shells (AA/ - = 2). It exists in 
clusters of any shape and exhausts most of the quadrupole strength. Though the resonance is 
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FIG. 2: Quadrupole photoabsorption in light deformed clusters Na^ , Na^ and Na^ 5 . The results 
are given as bars for every discrete state and as smooth strength functions obtained by folding 
with a Lorentzian of width A =0.25 eV. The weak infrared part of the strength (enclosed by the 
circles at 0.5-1.5 eV) is rescaled by the factor 2 • 10 2 . 
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FIG. 3: The electron level schemes for Na^ , Na^ and Na^~ 5 in the spherical limit (left) and at 
the equilibrium deformation (right). Occupied and unoccupied single-particle levels are drawn by 
solid and dash lines, respectively. The Fermi (HOMO) level is marked by index F. Arrows depict 
the possible low-energy hole-electron E2fj, transitions. Double arrows represent transitions of the 
scissor nature between the members the deformation multiplet. 

energetically close to the dipole Mie plasmon, it can be discriminated by means of angular- 
resolved electron energy-loss spectroscopy (AR-EELS) at electron scattering angles ~ 6° 

In the present study, we are not so much interested in the quadrupole plasmon but in 
non-collective infrared quadrupole modes which can deliver important information about 
the electron single-particle spectrum near the Fermi (=HOMO) level. The infrared spectra 
are associated with the low-energy AJ\f = branch created by E2 transitions inside the 
valence shell. Being of the AJV = origin, most of these spectra can exist only in clusters 
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with partly occupied valence shell, i.e. in deformed clusters. The deformation splits the 
infrared quadrupole mode into branches Xfi =20, 21 and 22. In Fig. 1, the infrared modes 
reside at 0.5-1.5 eV. As compared with the quadrupole plasmon, the infrared spectrum is 
very dilute. It is represented only by a few well separated levels. This prevents mixing of 
leh configurations by the residual interaction and creation of collective states. The infrared 
quadrupole modes persist to keep their leh nature. As is seen from Fig. 1, they have very 
weak quadrupole strength in the photoabsorption. But they may be accessible in two-photon 
processes. 

Fig. El shows single-particle levels and E^/z-transitions inside the valence shells in Nat, 
Na^, and Na^ 5 . The levels can be characterized by Nilsson-Clemenger quantum numbers 16| 
v = [Mn z N\. The leh pairs corresponding to the transitions read {[101] — [110]}2i in Na^ , 
{[220] - [200]} 20 , {[220] - [211]} 21 , {[220] - [202]} 22 in Na+, and {[220] - [200]} 20 , {[211] - 
[202]} 21 , {[211]-[200]} 2 i, {[220]-[202]} 22 inNa+. Following this scheme, the infrared modes 
in Fig. 121 can be unambiguously identified as particular leh configurations. The modes have 
different photoabsorption strengths depending on value of their E2-transition matrix element 
(le/z|r 2 Y" 2At |0). The RPA calculations which allow in principle any composition of of states 
confirm that the infrared modes are indeed almost pure leh states. In every mode the 
dominant leh component typically attains 99 — 100%. 

The infrared modes provide valuable information about cluster's properties. As is seen 
from most of the transitions connect the levels arising due to deformation splitting. The 
corresponding infrared modes are determined by the deformation and vanish at the spherical 
shape. Thus they deliver information on the deformation splitting of the electron levels in 



the HOMO-LUMO region. Besides, 
of occupied electron states (see e.g. 



aeing combined with photoemission data for the spectra 



24j). the infrared leh modes immediately yield energies 



of unoccupied electron states. As a result, one can get the complete electron spectrum in 
the HOMO-LUMO region. This spectrum is sensitive to different cluster features (interplay 
of equilibrium and isomer cluster's shapes, ionic structure, correlations, temperature effects, 
etc.) and so can serve as an effective tool for investigation of these features. 

Finally, it is worth noting that in deformed systems the electric and magnetic modes 
with the same projection /i and space parity tt are mixed. The mixture of electric E21 and 
magnetic orbital Mil excitations is especially interesting as it provides access to the orbital 
Ml scissors mode P, J, The properties of this mode are sketched in the next subsection. 
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FIG. 4: Macroscopic view of scissors mode : rigid rotation \2nS (left), and rotation within a rigid 
surface [3] (right). 



D. Scissors mode 



The scissors mode (SM) is a general collective flow already predicted or found in different 
finite quantum systems (atomic nuclei 



251 l26j | , quantum dots |27[ and ultra-cold superfluid 



gas of fermionic atoms [28j). Besides, this mode is used as an indicator of Bose-Einstein 



condensate in the dilute gas of trapped Bose atoms 



predicted 



29 



301 ] . In atomic clusters, the SM was 



but, in spite of some attempts j3l|, not yet observed. Obviously, the SM 
demonstrates a universal character. It is pertinent to any finite systems with two features 
in common: broken spherical symmetry (deformation) and a two-component nature. In the 
systems listed above, these components are neutrons and protons in nuclei, valence electrons 
and ions in atomic clusters, electrons and surrounding media in quantum dots, atoms and 
the trap in dilute Fermi and Bose gases. The SM features were recently reviewed in j6j. 



1. Macroscopic view 

The macroscopic collective nature of SM can be illustrated by the geometrical model first 
proposed for atomic nuclei [23. For clusters, this model exhibits SM as scissors-like small- 
amplitude oscillations of valence electrons versus the ions, both assumed to form distinct 
spheroids (see left part of Fig. 0]). Hence the name scissors mode. Obviously, such a mode 
can exist only in deformed systems. 

Following the alternative view of [3| (right part of Fig. Hj), the displacement field of the 



mode is a sum of the rigid small-amplitude rotation and irrotational quadrupole term (the 
latter provides vanishing velocity of electrons at the surface): 

u(r) = Q xr + 5 2 (l + 5 2 /3)- 1 V(yz) (1) 

where Q is the angular velocity and 5 2 = 3/2 — Rj_)/(R% + 2.Rj_) is the parameter 
of quadrupole deformation exp ressed via the semi-axes of the ellipsoidal system. Both 
macroscopic treatments, 0, include the rigid rotation of valence electrons versus the 
ions with the restoring force originating from the Coulomb interaction between the electrons 
and ions. But the model 3] has the additional irrotational quadrupole term. Its physical 
sense becomes clear if we take into account that the SM, like the quadrupole E2 modes, 
is separated into low-energy (AAf = 0) and high-energy (A7V = 2) branches. The first 
rotational term in (|T} is responsible for the AJ\f = branch while the second irrotational 
term generates the AJV = 2 branch. Thus the branches have quite different physical nature. 
This feature is exploited as a signature of the Bose-Einstein condensate in dilute gas of 



trapped Bose atoms 



29 



301 ] . Unlike the normal phase which exhibits both SM branches, 
the condensate phase supports only the irrotational flow. So, a vanishing AAf = branch 
serves as a fingerprint of the condensate. 

The rotational low-energy branch AJV = corresponds to a large extend to the scissors 
scenario exhibited in left part of Fig. 0] and thus is usually treated as a true SM. In atomic 
clusters, just this branch is responsible for van Vleck paramagnetism and other effects of 
orbital magnetism. In what follows, we will concentrate on this part of SM. 

In axial clusters, the low-energy SM is generated by the orbital momentum fields L x and 
L y perpendicular to the symmetry axis z. The mode is represented by the states (A 71 " = 1 + > 
where A is the eigenvalue of L, and 7r is the space parity. Energy and magnetic strength of 



the mode are estimated as 



, 3. 4, 61 



" = ^-K ll3 h eV, (2) 

s 

B(M1) = 4(1+ | L x | 0) 2 fi 2 b = 2 -Npu^ b ~ N^ 3 5 2 £ (3) 

where N e is the number of valence electrons, r s the Wigner-Seitz radius (in A), and fib is the 
Bohr magneton. We use here atomic units m e — H — c — 1. The value B(M1) stands for 
summed strength of the degenerated x- and y-branches. The z-branch vanishes for symmetry 
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reasons. It is worth noting that B (Ml transition strength does not depend on r s and so 
is the same for different metals. Both energy and magnetic strength are proportional to 
the deformation parameter 5 2 and thus vanish in spherical systems. So, the SM exists only 
in system with a broken spherical symmetry. The larger the deformation, the stronger the 
low-energy scissors mode. 

It worth noting that we assume here global deformation of the system. At the same time, 
spherical symmetry can be broken locally while keeping the global spherical shape. This 
takes place, e.g., in spherical clusters where the ionic background as such destroys locally 
the spherical symmetry and gives rise to some Ml (though rather weak) orbital strength 
even with zero global deformation [32]. 

The irrotational AJ\f = 2 branch of SM takes place in both systems, spherical and 
deformed. Its energy maps the energy of the quadrupole plasmon. In fact the irrotational 
SM is a part of this plasmon. 



2. Microscopic view 

The microscopic treatment of the SM is based on the shell structure of axially deformed 
mean field Q i Tfie angular _ o rt fio gonal t o tfie Symm e try aMs, 4 and 
promote low-energy AjV = transitions inside the valence shell and high-energy AM = 2 
transitions across two shells. 

It is instructive to expand the wave functions of the single-electron states in terms of the 
spherical basis (nLA) 

*,=[A,n,A] = Yl < L RnL(r)Y LA (n). (4) 
nL 

This allows to evaluate the single-particle orbital Ml transition amplitude between hole 
[y = h) and particle [y = p) states: 

(%\L x \* h ) oc 8 w 5 ApA ^J2<L< L VL(L+l)-A h (A h ±l). (5) 

nL 

Eq. © shows that the scissors mode is generated by A p = ± 1 transitions between the 
components of the same spherical (nL)-level. In spherical systems (nLA)-states belonging 
to the level (nL) are degenerate while in deformed systems they exhibit the deformation 
splitting and so may be connected by Ml transitions with non-zero excitation energies. This 
is the microscopic origin of the scissors mode. The energy scale of the scissors mode is 
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FIG. 5: E21 and scissors Ml photoabsorption for infrared modes in light deformed clusters. 

determined by the deformation energy splitting and so is rather small. This explains the 
predominantly low-energy (AAf =0) character of the true SM. 

There is an intimate connection between the scissors and quadrupole E21 modes in de- 
formed clusters. If fact, both SM and E21 are parts of a general motion characterized by the 
states \A n = 1 + >. To illustrate this point, let's again consider Eq. (JHJ). As was mentioned 
above, the SM operator connects only components from one and the same basis level nL. 
This is because the operators L x and L y have no r-dependent part and so, due to orthog- 
onality of the radial wave functions i? n ^(r), cannot connect the components with different 
nL. But the latter can be done by the quadrupole operator r 2 F 2 i- In this sense the SM 
operator is more selective than E21, though both operators generate transitions of the same 
multipolarity. The states \A n = 1 + > involve both SM and E21 modes and respond to both 
Mil and E21 external fields. The states are treated as magnetic or electric, depending on 
each of the two responses dominates. 

The dual nature of the states [A 77 = 1 + > is demonstrated in Fig. |S| The one-to-one 
correspondence between the E21 and SM peaks takes place. These peaks represent leh states 
{[101] - [110]} 21 in Na+, {[220] - [211]} 21 in Na+, and {[211] - [200]} 21 , {[211] - [202]} 21 
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in NaJ" 5 . The first free states originate from the deformation splitting of spherical levels and 
thus give rise to SM. In principle, these states carry both E21 and SM flows and the ratio 
between two contributions depends on the actual state structure. The plot for Na+ 5 gives an 
instructive example. It is seen that the lower peak determined by [211] — ► [202] transition 
between the members of the deformation multiplet exhibits an appreciable magnetic dipole 
strength and thus should carry a large SM fraction. Just such deformation-induced infrared 
states provide access to SM. Instead, the higher peak is determined by a [211] — ► [200] 
transition which takes place even in spherical case. The deformation is not crucial here. 
This state favors the E21 field and so can be treated as an ordinary quadrupole mode. 



3. Effects of orbital magnetism 



Static orbital magnetism in clusters was widely explored during last decades. The stud- 



ies were mainly based on the Landau theory of the atomic magnetism 



A variety of 

n - 



tibility 



issues was covered: giant dia- and paramagnetism in weak magnetic fields |34l kluj . size 
and temperature effects 36,|s7|, manifestation of quantum supershells in magnetic suscep- 



influence of cluster shape (both axial and triaxial) on t 



re magnetic properties 
38[ , orbital magnetism 



3l. l38l| . anisotropy of mag netic susceptibility in deformed clusters [3, 
of supported clusters |3jj], etc.. 

In the present study, we will consider some principle points of the orbital magnetism, 
connected with the scissors mode: the decisive role of SM in van Vleck paramagnetism 
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and the related effect of dia-para anisotropy in magnetic susceptibility of particular 
light clusters p . The latter effect displays a peculiarity of small clusters to exhibit strong 
variations in their properties with a cluster size. For example, RPA calculations jj, (| show 
that the SM energies and B(M1) strengths, though mainly scaling with the deformation 82 
and the electron number N e according to the trends (J2J) and (JHJ), demonstrate, nevertheless, 
strong fluctuations. This can affect the magnetic susceptibility and lead, in particular cases, 
to dia-para anisotropy. 

At this point, we introduce a few equations to make the discussion more precise. The 
interaction of cluster valence electrons with a uniform magnetic field B k applied along the 
coordinate axis k is 



H, 



int 



\i h B k L k + -fi b B kPk 



(6) 
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where k = x, y, z is the coordinate index, L k = Y^a=i * s ^ ne ^-th projection of the 
angular moment operator (the sum runs over all valence electrons), p 2 = J2a=i( x a + Va)i 



p 2 x = ^a=i(^o + z a)' an d p\ — J2a=i( x a + z a)- We neglected in © electron spins since for 
clusters considered below (axial sodium clusters with even N e and completely filled Fermi 
level) their contribution to the magnetic susceptibility is expected to be small. 

If the magnetic field is weak, one can use the perturbation theory. Then, up to the second 
order to B k , the induced change of the ground state energy is 

<j\L k \0 > | 2 



u 



int 



fM b B k <0\L k \0>-fM 2 b B 2 k J2 
1 



+ ^ 2 b B 2 k N e p 2 (7) 

where Uj is the energy of the excited state \j > and pi is the average value of p\. The negative 
second and positive third terms in (J2J are responsible for the temperature independent van 
Vleck paramagnetism and Langevin diamagnetism, respectively. 

The first (linear) term in (JJJ) dominates in the systems with a partly filled Fermi level 
since in this case there is only incomplete mutual compensation of the contributions of the 
valence electrons with different orbital projections A. In particular, this term results in a 
strong diamagnetic moments p = p b \M in axial clusters with odd N e 38]. In magic spherical 
clusters, where the Fermi level is fully occupied, both the first linear and second quadratic 
terms are zero and these clusters are again diamagnetic [Ssj] . This diamagnetism is called 
giant (3^ since, due to p\ ^> a 2 , (a is the Bohr radius), it is much stronger than the atomic 
one. 

We will consider axially (z-symmetric) deformed clusters with even N e and fully occupied 
Fermi level. In this case, the linear term in (JJJ) is zero but there remains the van Vleck term 
for k = x,y. The orbital magnetic susceptibility Xk — —d 2 oj™ t /dB 2 is then the sum of 
Langevin diamagnetic and van Vleck paramagnetic terms: 

Xk = Xt + X P r, (8) 

where 



X f a = -p 2 b N eP 2 > = -p 2 b Q§, (9) 

x p r = m E I<J|L " |0>|2 = ^ 
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FIG. 6: Normalized diamagnetic, paramagnetic and summed moments fi = x&k (Bk = AT) in 
axial deformed clusters Na^, Na^ 5 , Na^, Na^~ 7 , Na^" 5 , Na^" 5 and Na^~ 19 . 



having denoted by 

the cranking moment of inertia and by 



©f = N eP l 



<j\L k \Q > 



UJ; 



(11) 



(12) 



its rigid moment of inertia. 

Note that for k = x,y the operator entering in the matrix element in (fTUj) is exactly the 
scissors generator. This makes evident that just the low-energy SM mainly contributes to 
Xx% a - Following our calculations this contribution achieves 85 — 100%. So, just the SM 
determines the van Vleck paramagnetism. 

In the schematic model [3|, the moment of inertia comes out as the rigid-body value, so 
that 9 X)V = 6^ y and — ~xt,yy i- e - a complete compensation of dia- and paramagnetic 
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FIG. 7: Macroscopic image of twist mode. 



terms in Xx,y takes place. Due to axial symmetry, one also has x p z ara = 0- The total 
susceptibility becomes, therefore, strictly anisotropic j^] 

X* = Xy = 0, Xz = xf a , (13) 

varying from zero to diamagnetic values. 

However, strong shell effects in some light clusters may alter appreciably the above result. 
This is illustrated for Na^~ 7 in Fig. El Because of very low excitation energy of the SM in this 
cluster (4 , |(| , the paramagnetic susceptibility is so much enhanced that it cannot be balanced 
by the diamagnetic term. So, Na^~ 7 gains the remarkable property to be paramagnetic in x,y- 
directions and diamagnetic in z-direction. The cluster Na^ also hints this property. Though 
dia- and paramagnetism are weak as compared with other forms of magnetism, the dia-para 
anisotropy described above is big enough to be measured experimentally. This effect is only 
one of the examples of diverse unusual magnetic properties of small clusters (see j^] for a 
recent review). Altogether, these properties may provide an interesting perspective for new 
technologies. 



E. Twist mode 



The SM is the dominant orbital magnetic mode in deformed clusters. But in spheri- 
cal clusters the SM vanishes and the next orbital magnetic mode, twist M2, becomes the 

h 

strongest one pj. 

This twist was first proposed as an quadrupole torsional vibrational flow of an elastic globe 
j^j]. The mode is generated by the operator T = e~ tazLz = e au '^ with the velocity field 
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u = (yz, —xz, 0) and small angle a |40l l41|. The macroscopic image of twist is presented in 
Fig. The mode is treated as small-amplitude rotation-like oscillations of upper hemisphere 
against the lower hemisphere. The rotational angle of horizontal layers is proportional to z 
(projection to the axis of rotation). So, the flow vanishes at the equator and poles of the 
system. 

The restoring force of the twist is determined by the quadrupole distortions of the Fermi 
surface in the momentum space. So, the twist represents transverse magnetic quadrupole 
oscillations of an elastic medium, generated by variation of the kinetic-energy density. The 
twist is a general feature of any 3-dimensional finite Fermi system which exhibits an elastic 
behavior. The twist mode is well known in atomic nuclei and still remains an actual topic 



of both theoretical and experimental studies [41]- 



44j . The interest on twist is motivated 



by the fact that this mode is a remarkable example of an elastic vortex motion. The twist 
was also predicted but not yet observed in atomic clusters 0, and trapped dilute gas of 
Fermi atoms j^J. 

Unlike the Ml scissor mode which exists only in deformed systems, the twist takes place 
in Fermi systems of any shape, spherical and deformed. As the magnetic quadrupole mode, 
the twist is represented by \A* = 2~ > states of unnatural space parity. Twist energy and 
M2 strength in atomic clusters are estimated as P] 



17eV A 2 r; 2 N~ 1/3 , B(M2) = 0.52r*N*f4 (14) 



where B(M2) is the probability of M2 transition from the ground state to 2~ twist state. It 
is easy to see that the generating operator of the twist, zL z oc r(Y w L z ), coincides with the 
orbital component fj, — of M2 transition operator |46( 

F(M2,fi) = f i b VWr[g s {Y 1 s} 2fl + ^{Y.L}^} (15) 

which is the sum of spin and orbital contributions with g s = 2 and g\ = 1. So, it is natural 
to consider the twist as a part of the orbital M2 motion. 

Let's consider now the microscopic features of the twist. Figure |S] exhibits results of the 
microscopic calculations |5j for the orbital M2 strength B(M2) = |(2-|F(M2)|0)| 2 in light, 
medium and heavy spherical sodium clusters. Only the orbital part of the M2 transition 
operator is used. It is seen that the dominant twist strength is concentrated in a single leh 
peak with the lowest energy. This peculiarity persists independent of N e , though the degree 
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FIG. 8: The distribution of M2-strength in spherical Na clusters of a different size. 
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FIG. 9: Single-particle level schemes for the mean held of valence electrons and M2 lh-le transitions 
in Na^ and Na^ 3 . The levels are labeled by {nl} where n is the number of radial nodes in the 
single-electron function and I is the orbital moment. Transitions corresponding to the dominant 
twist peak are depicted by bold arrows. Contributions of these transitions to the complete B{M2) 
strength are given in %. 

of concentration decreases with N e . The twist peak exhausts 100%, 80%, and 60% of the 
total M2 strength in Nag , Nag 3 , and A^a^ 315 , respectively. As was shown in this peak 
corresponds to the nodeless twist mode plotted in Fig. [7| 
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The microscopic nature of the twist peak is exemplified in Fig. It represents the 
n, I — > n, I + 1 transition between single-particle levels with the node number n — 1 and 
maximal orbital moments /. The levels belong to the last occupied and first empty electron 
quantum shells. As a result, the twist can in principle deliver information about single- 
particle levels with maximal orbital momenta near the Fermi surface. 

Calculations P show that the twist strength dominates over the spin strength already 
in clusters with N e = 40. Starting with N e = 92, the orbital contribution becomes over- 
whelming and demonstrates for N e = 440 a huge value of 2 • 10 5 /i^A 2 eV. As both spin 
and orbital long-wave Ml responses are forbidden in spherical clusters (their Ml-transition 
matrix elements are proportional to the radial integral f Rn^^R^i^dr = 5 ni i 1>n2 i 2 and 
thus vanish in the non-diagonal case subject to ortho-normalization requirement), the twist 
becomes the strongest multipole magnetic mode already for clusters of medium size. This 
emphasizes its fundamental character. 

The search for the twist in atomic clusters is a demanding problem. This mode cannot 
be observed in photoabsorption (only dipole) or two-photon processes (only natural parity). 
Exploration of twist in atomic nuclei hints that one might try to search this mode in clusters 
by means of angular resolved electron-energy-loss-spectroscopy (EELS) 
twist contributes to the transverse form-factor, the back angles of scattered electrons are 
most suitable. 



23, 



47, m. Since 



III. EXPERIMENTAL ACCESS TO NON-DIPOLE MODES 



A. General view 



As compared with other non-dipole modes, infrared quadrupole states deserve special 
attention. There are at least two samples, i) free light deformed clusters and ii) oriented 
silver rods embedded in a matrix, where these modes can be studied experimentally in 
two-photon processes. The light deformed clusters are especially interesting j?J since their 
quadrupole modes provide access to the single electron spectra. In turn, the spectra can 
serve as a sensitive indicator of several cluster properties. It is known that just small 
clusters demonstrate properties unusual for the bulk and, in this sense, are attractive for 
both fundamental physics and applications. 
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Two-photon processes (TPP) are widely used in atomic and molecular physics (for a 
comprehensive discussion see |52|). Since atomic clusters are similar to molecules, it would 
be natural to implement the same reactions to clusters. However, applications of TPP 
to clusters are still very scarce. This can be partly explained by our poor knowledge on 
non-dipole low-energy electron modes in clusters, and partly by peculiarities of clusters. 

In this section we discuss applicability of some traditional TPP methods to clusters 
with particular emphasis on infrared quadrupole modes. In this connection, some cluster 
properties essential for TPP should be closely scrutinized. Namely: 

i) The dominant decay channel of cluster levels is usually not radiative. The collective 
states, like the dipole plasmon, mainly decay through Landau fragmentation (dissipation of 
the collective motion through surrounding leh excitations). For non-collective states, the 
electron-electron collisions and electron-ion coupling are most essential. 

ii) Collective electron modes in clusters have very short lifetimes. The typical example is 
the dipole plasmon which represents in medium large clusters with N e ~ 50 — 1000 a broad 
( ~ 0.5-1 eV) mode with a life time of 10-20 fs. This is due to a coupling of the collec- 
tive strength with energetically close leh states whose availability increases with increasing 



cluster size 



49, 



50L_The mechanism is similar to Landau damping of the plasmon in the 
. Its analogue in finite systems is called therefore Landau fragmen- 



bulk electron gas 

tation. The short lifetimes hamper application of adiabatic TPP methods as they would 
require more intense laser irradiation. The intense irradiation results in strong dynamical 
Stark shifts and undesirable population of high-energy states by second harmonic of the laser 
pulses. These effects can be detrimental for adiabatic TPP especially if the cluster spec- 
trum is not sufficiently dilute. Non- collective electron modes (as low-lying quadrupole leh 
ones) have much longer life times that can approach ps for low temperature and hundreds 
fs for the room temperature 12, Q|. Unlike molecules, these modes in metal clusters are 
mainly driven by the mean field similar to the nuclear one. This allows their unambiguous 
treatment at least in light clusters with 5-20 atoms. 

iii) Clusters exhibit a strong shape isomerism. Even beams with size selected clusters 
usually represent a statistical mix of samples with a multitude of different shapes. This 
blurs the measured low-energy electron spectra. 

The properties i)— iii) are absent in atoms but at least some of them are common for 
molecules. Thus one may hope that methods of molecular spectroscopy can be applied 
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FIG. 10: Two-photon process: scheme of population of infrared quadrupole modes \[i =20, 21 
and 22 via the dipole states with \[i = 10 (left) and 11 (right). 



to clusters. It is worth noting that, as compared with molecules, clusters provide new 
possibilities for investigation of electron modes. The electron spectra in clusters are very 
diverse (from the collective plasmons to the pure leh configurations) and, at the same time, 
are easily classified and treated in terms of electronic quantum shells. 

Light free clusters with 5-20 atoms have dilute spectra, a feature which allows to avoid 
some of the troubles listed above. Landau fragmentation in such clusters is very weak, the 
infrared levels are not so broad and their lifetimes are much longer. Most of light clusters 
have larger isomer energies such that they have one clearly preferred equilibrium shape and 
many of the energetically close isomeric shapes are similar to the ground state one 22j, thus 
minimizing blurring of the low-energy spectra. Beams of size-selected singly- charged light 
clusters are readily available. The cluster temperature ~ 100-200 K could be optimal for 
our aims. Then the thermal broadening of electron levels is small. 

Oriented si.ver rods embedded into glass Q Q represent another interesting samp.e to 
be studied in TPP. These clusters are known as commercial polarizers. The polarization is 
provided by the property of the A/i = 10 branch of the dipole plasmon to absorb only the 
light aligned along the symmetry axis of the cluster. Because of the big axis ratio, the rods 
represent an example of atomic clusters with an extreme deformation. These heavy clusters 
can be studied in Raman scattering jloT ]. 

Following the general TPP scheme presented in Fig. QUJ the target quadrupole level is 
populated via an intermediate dipole state by two (absorption and emission) dipole transi- 
tions. Both the collective dipole plasmon and non-collective leh configurations can serve as 
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FIG. 11: TPP population transfer from the ground state |0 + > to the quadrupole state |2 + > 
via the leh dipole state 1 1 > in Na^. The detuning from the dipole state is A ~ 0.14 eV. 

intermediate states. The TPP via the dipole plasmon, with its broad and smoothed structure 
and ultra-fast decay through Landau damping, reminds the population transfer via contin- 
uum. The TPP via non-collective leh dipole states is more transparent and represents the 
typical population transfer process in tree-bound-state A-systems. 

As is shown in Fig. EH there are possible TPP paths via dipole states with A/i = 10 and 
11. These states could be, for instance, branches of the dipole plasmon, separated by the 
deformation splitting. In the first case, the population of A/i =22 mode is forbidden and 
only 20 and 21 modes can be targeted. This provides some selection when carrying out TPP 
experiments. 

For a crude analysis of TPP in clusters, one may use the rotating wave approximation 
(RWA) that allows to discard the fast laser frequencies and treat the problem with the 
much lower Rabi frequencies only. Fig. ^2 for Na^ shows a typical electron A-system 
and hence the scale of the laser frequencies corresponding to photoabsorption (pump) and 
photoemission (Stokes) processes. Following our calculations, even for laser pulses with 
intensity / ~ 10 l0 W / cm 2 , Rabi frequencies in light deformed sodium clusters lie in the 
interval 10 _1 — 10~ 2 eV. Obviously, they are much less than the carrier (laser) frequencies 
and thus fulfill the main RWA requirement. Besides, because of the pretty long life-time of 
the target non-collective quadrupole states, one may operate with long laser pulses covering 
many Rabi cycles. Hence we match the second RWA requirement. 

In the next subsections, we will consider some widespread two-photon processes (Raman 
Scattering (RS), Stimulated Emission Pumping (SEP), and Stimulated Raman Adiabatic 
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Passage (STIRAP)) and estimate their ability to observe infrared quadrupole modes. 



B. Raman scattering 

RS is one of the simplest two-photon processes where a dipole laser-induced transition 
to an intermediate electronic level (real or virtual) is followed by the dipole spontaneous 
fluorescence to low-energy levels. 

RS measurements of electron infrared modes in clusters are yet very rare. For instance, the 
quadrupole infrared modes were observed in heavy silver clusters embedded into amorphous 
silica films [31]. The dipole plasmon was used as a resonant intermediate state. In the 
recent experiment of this group with another sample (oriented silver rods embedded into 
glass), the infrared electron quadrupole mode, probably A/i = 20, was detected [lOj. These 
observations show that the dipole plasmon can be used as intermediate state for TPP in spite 
of its extremely short lifetime. This is the encouraging message for TPP implementation to 
clusters. 

RS generally assumes that the coupling between the intermediate dipole and target 
quadrupole states, from one side, and the coupling between the intermediate and ground 
states, from the other side, are of the same order of magnitude. Only then the final state is 
successfully populated. This means that the absorption and emission dipole matrix elements 
should be of the same scale. Experiments of E. Duval Q, [yj show that this is indeed the 
case for heavy clusters. In light deformed clusters the situation is more complicated. Our 
calculations |7[ for Na^", Na^ and Na^" 5 revealed that the absorption and emission dipole 
matrix elements lie mainly in the interval 2-10 eao (atomic units) and thus are basically sim- 
ilar. This means that infrared quadrupole modes in light deformed clusters have a chance 
to be observed in the resonant RS running via the dipole plasmon. However, we have also 
found a few exceptions when the Stokes couplings are strongly suppressed due to destruc- 
tive interference effects (see discussion in j?J). The particular example is oblate Na^ where 
only the path via the A/i = 11 branch of dipole plasmon results in a strong TPP while the 
alternative path via the A/i = 10 branch is about completely suppressed. Such considerable 
variations for different paths are not surprising for small clusters where electronic spectra 
and transition rates are very specific and can appreciably change from one sample to an- 
other. In clusters with suppressed dipole couplings, the RS is not effective. Besides, the RS 
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has a general drawback that its photoemission step is fully dictated by competition between 
different decay channels of the intermediate state. As a result, the spontaneous fluorescence 
to the target state and hence the RS efficiency is typically rather low. To overcome this 
trouble, we should consider methods with stimulated emission into the target level. 



C. Stimulated emission pumping 

The stimulated emission pumping (SEP) enjoys widespread application for atoms and 
molecules and seems to be very promising for population of infrared quadrupole states in 
clusters. Unlike RS, this method exploits two lasers, pump and Stokes (or dump) [53, 



The pump pulse is responsible for the first photoabsorption step. The Stokes pulse couples 
the intermediate state to the target one. The pulses can be simultaneous or sequential. If 
the difference between the pump and Stokes frequencies is in resonance with the frequency 
of the target state, then the Stokes radiation stimulates emission to this state. 

One should distinguish SEP with incoherent and coherent irradiation. If incoherent light 
sources are used, then, even in the best case of the saturated process and coincident pump 
and Stokes pulses, the maximal SEP efficiency achieves only one-third of the complete 
population 56] . In practice, the transfer efficiency does not exceed 10%. However, even this 
value is enough for many spectroscopic studies. 

A much better population can be achieved for coherent laser irradiation. If the pump 
and Stokes pulses coincide, then under certain requirements (given in the next subsection) 
we can even gain the fully adiabatic process (with still not complete population transfer as 
compared with STIRAP j^). In this subsection, we consider the coherent SEP with an 
appreciable detuning (often called as off-resonant stimulated Raman). This process promises 
a simple, robust and efficient population scheme. 

We present below our TDLDA results for SEP in Na^. The calculations explore the 
time evolution of the initial single-electron states. All the possible channels including pho- 
toionization (through absorbing boundary conditions) are taken into account. The ions 
are considered in soft jellium approximation at fixed deformation. The details of the 
method are described elsewhere [id Q]. The parameters of the Stokes and pump pulses 
are given the panel (c). Their frequencies maintain the two-photon resonance with the tar- 
get state: u p — oj s = co> 2 +o and are detuned somewhat below the intermediate dipole state, 
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FIG. 12: Off-resonant SEP population transfer from the ground state to the quadrupole state 
|2 + > via the leh dipole state | X > in Na^. The left plots exhibit (in arbitrary units) strengths 
of quadrupole (a) and dipole (c) modes as a function of their energy. The right plots depict (in 
atomic units) time evolutions of the overwhelming quadrupole mode with w 2 +=0.80 eV (b) and 
dipole oscillations (d). Specifically, the oscillations of the quadrupole and dipole moments are 
presented. The parameters of the Stokes (s) and and pump (p) pulses (frequency u, duration T 
and intensity I) are given in the panel (c). The pulse profile is the squared cosine. 

A = {uji- — uj q +) — u p = (loi- — u 2 +) — uj s = 0.14 eV, so as to minimize its population. It is 
worth noting that an appreciable detuning is crucial as the dipole couplings in clusters are 
rather strong. This is the only way to suppress undesirable population of other dipole and 
quadrupole states and thus to enhance the population of the target state. 

In Fig. the TPP via the isolated non-collective dipole level at 1.38 eV is considered 
(the actual A-system is presented in Fig. HT| . The quadrupole (A = 2) and dipole (A = 1) 
strengths <tea(^) = SjO l rA ^A/jO) 2 5(a;j— u) versus their excitation energy are depicted in the 
left panels. Panel (a) demonstrates the strong quadrupole excitation at 0.80 eV. Because of 
the considerable detuning, all other quadrupole and dipole modes are pretty small. The right 
panels (b) and (d) show that, while dipole oscillations survive only during the irradiation 
time ~ 200 fs, the quadrupole oscillations around the static values Q2 = 26.6 clq (ao is the 
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FIG. 13: Time evolution of dipole and quadrupole oscillations for different pulse shifts (sequences) 
in Najp The values of the shifts between the Stokes (s) and pump (p) pulses are depicted at the 
left plots. The pulse parameters are the same as in Fig. ^Jfor exception of I s = 3.3 • W 10 W/cm 2 . 

Bohr atomic unit) are very endurant. There is no any attenuation until 800 fs. Our study 
shows that the lifetime of this mode is even longer and achieves at least a few ps. Because of 
the dilute spectrum, decay of this mode via the electron-electron coupling is negligible. Only 
the coupling with the ions remains but it is not taken into account in these calculations. 
The electron-ion coupling should limit lifetime of the mode by a few hundred fs at a room 
temperature and by a few ps at a minimal temperature ~10 K. The latter case corresponds 
to our analysis. 

The calculations show that the endurant quadrupole mode is also produced in TPP 
running through the dipole plasmon. In this case, the mode is even several times stronger. 
The pulses with the frequencies up =2.04 eV and up =1.24 eV and the same intensities and 
durations as in Fig. H21 were used. The pulses were appreciably detuned from the dipole 
plasmon so as to weaken the competing channels (compare up =2.04 eV with the dipole 
plasmon in Na^ depicted in Fig. [TJ). It is remarkable that the dipole plasmon can be used as 
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the effective mediator in spite of its extremely short lifetime. The key point is the detuning 
that prevents the population of the plasmon and thus the leaking from it. Using the dipole 
plasmon as the intermediate state has the additional advantage of dealing with the well 
know structure rather than with still poorly known low-energy dipole states. 

Fig. 1121 demonstrates the dependence of the quadrupole signal on the sequence of the 
pump and Stokes pulses. It seen that the quadrupole mode is maximal at simultaneous 
pulses. This is the typical SEP case. The process is rather robust to changes of pulse 
characteristics. The calculations show that variations by several times in pulse durations 
and intensities still maintain the stable and distinctive SEP (though with the rescaled mode 
amplitude). Of course, any change in the pulse intensities should be compensated by the 
corresponding modification of the detuning. 

It is a non-trivial problem to detect the quadrupole mode populated in SEP. In principle, 
this could be done by using the probe laser with supplemented PES (photoelectron spectra) 
spectroscopy. The probe laser should be intense enough to result in measurable PES and its 
frequency should be large to map the single-electron states in the cluster and thus provide 
the complete PES [58]. In Na^, it could be u> pro be ~ 7 - 8 eV. The probe pulse should 
follow the pump one with a considerable delay so as to avoid overlapping with the latter and 
map only the endurant quadrupole mode. Due to the Raman coupling of the single-electron 
PES with the quadrupole oscillation, every PES should acquire left and right satellites with 
the relative frequencies coinciding with uj 2 +. Hence we get the frequency of the quadrupole 
state. The strategy is to keep the fixed pump frequency with a reasonable detuning bellow 
the intermediate dipole state and scan the Stokes frequency. When we get the resonance 
satellites up — u s = u>2+, then the endurant quadrupole mode has to be produced and 
detected. It may happen that the Stokes frequency maps a dipole state. In that case 
the endurant dipole oscillation are accompanied by a subsequent considerable increase of 
photoemission yield together with the appearance of dipole satellites. Such cases can be 
easily separated from the case of interest, population of the quadrupole mode, since the 
latter is not associated with a significant resonant increase in photoemission yield but only 
by the satellites signature alone. 

The SEP experiments usually need three lasers: pump, Stokes and probe. However, 
one may propose the simpler recipe with implementation of only pump and probe lasers. 
Such scheme assumes the combination of the direct two-photon excitation of the quadrupole 
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state and PES techniques. Specifically, the intense pump laser should scan until its two 
photons excite the quadrupole state. Then the delayed probe laser is applied in the manner 
described above to detect the population. The calculations show that this simple method is 
quite effective and robust. 



D. STIRAP 

The STIRAP method QSQ promises up to 100% population transfer from the initial 
to the target level. This high efficiency is achieved by the coherent adiabatic character of the 
process and the principle possibility to avoid any leaking from the intermediate levels. Like 
SEP, STIRAP also implements pump and Stokes lasers providing the couplings |0 > — 11 > 
and |1 > — 1 2 >, where |0 >, |1 > and |2 > are initial ground state, intermediate dipole 
state and target quadrupole states, respectively. However, STIRAP is much more involved. 
Its main requirements are: 

i) Two-resonance condition up — us = ^2 — ^0 with a possible detuning A = (u\ — ujq) — 
up = (u>i — UJ2) — ojs from the intermediate state frequency. 

ii) Counterintuitive sequence of partly overlapped pulses when the Stokes pulse comes 
before the pump one. 

iii) Adiabatic passage of the state vector from the initial to the target state. The adia- 
baticity condition (for the detuning A=0) has the form 

QAt » 1 (16) 

where £1 = v/H| + is the average of the pump and Stokes Rabi frequencies and At is the 
overlapping time of the pulses (duration of the evolution process). 

Under these requirements, one of the adiabatic time-dependent eigen-functions of the 
system becomes a superposition of the initial and target bare states only: \b (t) >= c (t)|0 > 
+C2(t) 1 2 >. This is so called dark state. It is reduced to |0 > at the beginning of the adiabatic 
evolution and to |2 > at the end. Hence, the system finally finds itself in the target state. 
The intermediate state |1 > is not involved in \bo(t) > and thus is not populated at any time. 
Thus any leaking in the population is avoided and the transfer is complete. The main point 
is to evolve the system adiabatically, keeping it all the time in the state \bo(t) >. STIRAP 
is well known in atomic and molecular spectroscopy. It is rather insensitive to precise pulse 
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characteristics. Both continuous and pulse lasers can be implemented. 

Let's now examine the STIRAP requirements for the case of infrared quadrupole modes 
in atomic clusters 

Two-resonance condition. The better fulfilled this condition, the less mixing of the dark 
adiabatic state \bo(t) > with the state |1 > and thus more robust the STIRAP. In practice, 
this condition is perturbed by the detrimental dynamical (time dependent) Stark shifts of 
the electron levels 

Too strong Stark shifts perturb the STIRAP balance and destroy the adiabatic process. The 
Rabi frequency reads |5j 



n= \AJE~ 2.20 • 10 8 |d[eao]| \//[-^] s- 1 , (18) 
ft v ceo V cm 

where d is the dipole coupling matrix element in atomic units, I is the laser intensity in 

W/cm 2 , c is the light speed, eo is the vacuum dielectric constant. It is seen that Q ~ v 7 ! 

and thus S ~ I. So, the smaller the laser intensity, the weaker the Stark shift. However, 

the intensity cannot be too small because of the adiabaticity condition (|16p. The detuning 

A also has the upper limit established by the Rabi frequency Q. 

Following our calculations, the coupling dipole matrix element in light clusters is d ~ 

1 — 5 eao- If to put the laser intensity I = 10 10 W/cm 2 , then Eq. ()18|) estimates the peak 

Rabi frequency as Q ~ 0.02 - 0.1 fs or hVL ~ 0.01 - 0.07 eV. These values are less than the 

detuning A=0.136 eV used for the SEP analysis in the previous subsection. Thus the SEP 

was obtained at the tail of the dipole coupling between the states of the A-system. Under 

the given detuning and laser intensity, the Stokes shift is S ~ 0.0004 - 0.01 eV which is still 

a bearable value. 

Counterintuitive order of pump and Stokes pulses. The order when Stokes pulse precedes 
the pump one is crucial for the complete population since only then the state vector can be 



fully projected into the dark adiabatic state \bo(t) > j54j, |55|. Stokes and pump pulses must 



overlap and the overlapping time At determines the duration of the adiabatic evolution. 

The At should not be longer than the lifetimes of the initial and target states while 
the lifetime of the intermediate state does not matter if this state is not populated. So, the 
ultra-short lifetime of the dipole plasmon (10-20 fs) is not the problem for the ideal STIRAP. 
Even under some leaking from the plasmon (and thus the modest coupling of the dark state 
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,(t) > with |1 >) STIRAP is still possible. This is partly confirmed by recent experiment 
I where STIRAP was induced via continuum. Even in this extreme case with ultra-fast 
decay channels and considerable leaking, a clear STIRAP population transfer was obtained 
(though with low efficiency ~ 6%). 

When STIRAP runs via an isolated non-collective dipole state, like that with uj = 0.80 
eV in Na^, then the situation is much simpler. Following the estimations for the lifetimes 
of leh state, the value of At should not exceed a few hundreds fs at room temperature and 
a few ps at low temperature. 

It should be emphasized that the maximal population transfer at the counterintuitive 
order of pump and Stokes pulses is the main signature of STIRAP. Just this feature allows 
finally to judge if we deal with STIRAP or other processes like SEP. 

Adiabatic passage. This requirement is crucial for STIRAP. Only adiabatic evolution can 
ensure the steady complete populations. The adiabaticity condition is given by Eq. (fTBj). If 
we put, following the above estimations, Q ~ 0.02 — 0.1/s and At ~ 300/s, then it is easy to 
see that the adiabaticity condition is fulfilled: QAt ~ 6 — 30 >> 1. However, the situation 
is not so simple since Eq. (116)1 does not take into account the detuning that effectively 
weakens the Rabi frequency. In practice, we should use an appreciable detuning to minimize 
the competing channels. But then Q is considerably decreased and the adiabatic condition 
is violated. This becomes a real hindrance when the electron spectrum is not dilute enough. 

It worth noting that the adiabaticity condition (|T6|) can be considerably relaxed. Indeed, 
it was derived for one intermediate level. But the realistic spectrum of the dipole plasmon 
consists of a broad structure consisting from a sequence of dipole levels (see Fig. 1). In 
this case, the STIRAP adiabatic condition should be revised. The realistic case is more 
complicated but, at the same time, opens new possibilities for STIRAP. In particular, one 
may be allowed to loose the adiabatic STIRAP condition and thus the requirements for the 
laser intensity. 

A general case of N intermediate states, each with its own coupling and detuning was 
studied in It was shown that the trapped adiabatic state |&o(^) > can be created 

only when the ratio between each pump coupling and the respective Stokes coupling is the 
same for all intermediate states. Following our calculations, this condition is unrealistic for 
atomic clusters. However, softer alternative adiabatic requirements can be formulated. In 
particular, in the general case of arbitrary couplings, one may tune the pump and Stokes 
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FIG. 14: Simultaneous resonant excitation of the dipole state and quadrupole plasmon in Na+ X by 
the first and second pump harmonics, respectively. The up and middle left plots exhibit strengths 
of the quadrupole and dipole modes as a function of their energy. The right plots depict time 
evolutions of these modes. 

lasers just below all intermediate states and thus form so called adiabatic-transfer state which 
also results in a high, although not complete, population of the target level. Unlike |6 (*) >, 
the adiabatic-transfer state can have admixture from the intermediate states during the 
evolution period At and so some population leaking is unavoidable. Nevertheless, we have 
here a solid adiabatic transfer with a high population of the target state. Just this method 
is implemented in our calculations while using the dipole plasmon as the intermediate TPP 
state. 

Altogether, the arguments and estimations presented above show that STIRAP is possible 
in clusters. However, it is not so trivial to find the optimal parameters of the process and 
get STIRAP in realistic TDLDA calculations. Serious problems arise in connection with 
dipole and quadrupole plasmons which reside in the energy region 2 - 4 eV and thus are 
always covered by multiphoton excitations of the pump and Stokes pulses. Hence there 
arises a considerable undesirable population leaking. Especially strong coupling with the 
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FIG. 15: Population transfer from the ground state to the quadrupole state |2 + > via the leh 
dipole state 1 1 > in Na^. The up and low panels represent the quadrupole and dipole strengths 
at different delays (shifts) of the pump pulse versus the Stokes one: 82 fs (counterintuitive STIRAP 
pulse sequence), fs (coinciding pulses), -82 fs (intuitive pulse sequence). The parameters of the 
Stokes (s) and and pump (p) pulses (frequency uj, duration T and intensity I) are given in the 
low-right panel. 

quadrupole plasmon takes place when the pump or Stokes frequencies are in resonance with 
some dipole state. The example of such process is exhibited in Fig. ITU where the resonant 
excitation of the dipole state at 1.38 eV leads to a huge population of the quadrupole 
plasmon. This means that we should avoid the resonance cases by imposing appreciable 
detuning. But, as was mentioned above, this will effectively weaken the Rabi frequencies. 
Then, to keep the requirement (JTBJ), one should use longer pulses. This is quite possible at 
low temperatures when the lifetimes of the states can reach a few ps. In this case, we even 
have the freedom to decrease the laser intensity with the aim to minimize the Stark shifts. 

The proper choice of the pulse parameters allows to produce a STIRAP-like process 
in clusters even at high laser intensity. A successful example is presented in Fig. 
A large detuning of A = 0.15 eV is here combined with a high pulse intensity of I ~ 
10 n W/cm 2 . This particular choice allows to minimize the impact of the quadrupole and 
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dipole plasmons and, at the same time, to keep the adiabaticity condition. Fig. [T3] shows 
that the counterintuitive pulse sequence provides the stronger strength of the quadrupole 
mode than the coinciding pulses. In terms of the population, this could be the unequivocal 
signature of the STIRAP. As is seen from the figure, the intuitive pulse sequence also results 
in the same quadrupole strength. It looks like this TPP embraces not only the typical 
STIRAP branch corresponding to the dark state | fc»o > but also the branch from the dark 
state | b + >. The latter takes place because of the large detuning below the intermediate 
dipole state. Folowing [6l]j], an appreciable detuning below or above the intermediate state 
results in the additional dark dressed state \b + > or |6_ >, respectively. 

Though the process in Fig. ^] looks like STIRAP, its STIRAP nature still should be 
proved. The point is that the strength exhibited in Fig. does not correspond to the 
population (squared amplitude of the target state) but rather to the coherence (product of 
amplitudes of initial and target states). Hence, the maxima of the strength at non-zero pulse 
shifts still do not mean that the similar maxima take place for the population of the target 
state. Our TDLDA model does not allow to estimate properly the population. However, this 
can be done within the Rotating Wave Approximation for 3- level A-system J^|. Following 
these estimations, we have found that the TPP in the example above has the maximal 
population at coinciding pulses and thus is not STIRAP. For getting the STIRAP we need 
more intense laser pulses and maybe another value of the detuning. One may get other, and 
more favorable, options in clusters by a more suitable choice of the process parameters, set 
to worked out. In any case, the example in Fig. EHhints the principle possibility of STIRAP 
in atomic clusters. The experimental detection of the STIRAP population can be done by 
the same methods as in SEP. 

The application of STIRAP to clusters seems to be generally feasible. At the same time, 
this process is more demanding to realize than SEP. It requires further studies within realistic 
TDLDA models to look for the most effective process configurations. STIRAP can hardly 
be more practicable than SEP in population and detection of infrared quadrupole modes 
in clusters. However, one may try to implement STIRAP to get more transfer efficiency in 
particular cases. Besides, STIRAP running via the dipole plasmon can be interesting for 
theory of adiabatic TPP as such. Indeed, the intermediate state like the dipole plasmon does 
not exist in atoms and molecules and so was not yet investigated. Perhaps, this intermediate 
state will offer new options for STIRAP and extend our knowledge about this intriguing TPP. 
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IV. CONCLUSION 



Non-dipole infrared electron modes in atomic clusters (electric quadrupole (E2), magnetic 
dipole scissors (Ml) and magnetic quadrupole twist (M2)) were analyzed theoretically using 
as tool the time-dependent local-density approximation (TDLDA). The collective and non- 
collective nature of these modes as well as their connection with the basic cluster properties 
(single electron spectra, orbital magnetism) were analyzed in detail. Free light deformed 
sodium clusters and embedded oriented silver rods were considered as the most promising 
samples. The possible routes of experimental observation of infrared quadrupole modes in 
two-photon processes (Raman scattering, stimulated emission pumping (SEP) and stimu- 
lated Raman adiabatic passage (STIRAP)) were scrutinized. Following our analysis, the 
combination of SEP and PES techniques provides probably the most simple and robust way 
to populate and detect the infrared quadrupole states in free light deformed clusters. As 
a complementary method, a technique using direct two-photon excitation was also consid- 
ered. The Raman scattering seems to be optimal for large clusters, especially embedded 
oriented silver rods. The calculations encourage applicability of STIRAP to clusters. How- 
ever, this intriguing but demanding method still needs further realistic TDLDA studies. 
Smaller clusters as, e.g., Na? are the most promising candidates. 

We hope that our study will stimulate investigation of non-dipole electron modes in 
clusters by means of TPP experimental methods. This would provide a deeper understanding 
both clusters properties and peculiarities of different TPP schemes. 
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